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ABSTRACT 


•^Measurements  of  the  acoustic  reflectivity  of  the  lower  surface  of 
sea  Ice  were  wade  at  several  location  In  the  channels  of  the  Canadian  Arctic 
Archipelago  from  1971  to  1973 • The  distance  between  the  projector  and  hydro< 
phone  was  typically  100  to  300  m;  the  gracing  angle  at  the  Ice-water  Inter- 
face ranged  up  to  the  frequency  Halts  were  200  Hz  and  20,000  Hz. 

The  measured  reflectivities^  which  are  presented  hereT^ show  large  excursions 
from  unity.  Reflectivities  of  -15  dB  are  not  uncommon  in  spite  of  the  uni- 
form flatness  of  the  ice-water  Interface  and  In  spite  of  the  shallowness  of 
the  grazing  angle.  However,  a calculation  based  on  Klrchhoff's  Integral 
Theorem  and  Involving  an  Integration  over  the  interface  indicates  that  the 
observed  undulation  in  the  Ice-water  interface  Is  sufficient  to  account  for 
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INTRODUCTION 


The  problem  of  understanding  and  predicting  the  propagation  of  sound 
in  the  ice-covered  seas  of  the  Arctic  requires  a knowledge  of  the  behaviour 
of  sound  waves  as  they  reflect  from  the  water-ice  boundary.  The  scattering 
or  reflection  of  sound  from  the  rough  wave-covered  surface  of  more  southern 
oceans  has  been  fairly  extensively  studied  by  (to  name  a few)  Liebermann 
(1948),  Rckart  (1953),  Marsh,  Schulkin  and  Kneale  (1961),  Marsh  (1963), 

Ferris  and  Kuperman  (1970).  Fortuin  (1970)  has  prepared  a survey  of  the 
pertinent  literature.  In  the  Arctic  seas,  however,  propagation  studies  have 
not  Included  the  detailed  study  of  the  reflection  of  sound  from  the  under 
surface  of  the  sea  ice. 

During  springtime  field  trips  in  1971,  1972, and  1973, short  range 
measurements  were  made  of  the  reflection  loss  of  acoustic  waves  under  the 
surface  of  the  sheet  of  ice  upon  which  the  camp  was  situated.  Although  the 
ice  surface,  and  therefore  presumably  the  under-ice  surface,  was  always 
quite  smooth,  the  reflection  losses  were  strongly  frequency  dependent  and 
as  great  as  25  dB  in  some  instances. 

This  report  describes  the  experiment  and  presents  the  measured 
reflectivities.  Also,  a model  based  on  Kirchhoff’s  Integral  Theorem  is 
developed  to  explain  the  large  variations  that  were  observed. 

LOCATIONS 

Measurements  were  made  in  April  1971  and  April  1972  in  Robeson 
Channel  and  in  May  1972  and  April  1973  in  Barrow  Strait  at  the  locations 
shown  in  Figure  1.  The  ice  on  which  the  camp  was  located  was,  in  all  cases, 
smooth,  first-year  ice  having  a thickness  of  1.8  m in  1971  and  1972  and 
2.4  m in  1973.  Figure  2 is  a photograph  showing  the  snow-covered  ice  surface 
in  the  neighbourhood  of  the  camp  in  Barrow  Strait.  In  1973,  with  the  aid 
of  Canadian  Forces  divers,  the  undulations  of  the  ice-water  interface  were 
mapped  with  a grid  spacing  of  1 meter  over  an  area  in  the  region  of  the  camp 
measuring  8 by  31  meters.  Figure  3 is  a photograph  of  the  underside  of  the 
ice  surface  and  Figure  4 is  an  isometric  diagram  of  the  mapped  surface. 

The  R.M.S.  variation  of  the  level  of  this  Interface  was  only  2 cm  and  the 
subjective  Impression  of  the  divers  was  that  the  interface  was  remarkably 
smooth. 


Figure  2.  Photograph  o£  an  ice-camp  Illustrating  the  smoothness  of  the  Ice 


A photograph  of  the  underside  of  the  ice  surface  illustrating  the 
smoothness  of  the  interface.  Projecting  into  the  water  and  giving 
a feeling  for  the  scale  of  the  picture  is  a bamboo  pole  and  flag. 
The  bubbles  at  the  surface  are  the  divers’  exhaust  air. 
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which  passed  a pre-set  number  of  cycles.  This  burst  was  amplified  and  trans- 
mitted to  the  projector.  The  resulting  acoustic  signal  was  detected  by  the 
hydrophone,  amplified,  filtered  to  Improve  the  signal  to  noise  ratio,  and 
displayed  on  the  oscilloscope.  This  sequence  was  repeated  every  few  seconds. 
The  signal  was  recorded  by  photographing  the  oscilloscope  face  with  a Polaroid 
camera.  This  technique  Is  similar  to  that  used  by  Ferris  and  Kuperman  (1970) 
and  Medwin  (1966). 


Figure  6.  Block  diagram  of  the  electronics. 

The  example  photograph  in  Figure  7 shows  the  time  separation  of  the 
signals  from  the  two  paths.  This  time  separation  sets  a lower  limit  to  the 
frequency  that  can  be  used  in  the  experiment  since  the  period  for  one  cycle 
must  be  substantially  less  than  the  time  delay  In  order  to  obtain  a reliable 
measurement.  In  order  to  make  measurements  at  lower  frequencies,  the  separa- 
tion of  the  hydrophone  and  projector  must  be  Increased  to  Increase  the  time 
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delay.  This  Increased  separation  unfortunately  reduces  the  signal  to  ambient 
noise  ratio  because  of  the  limited  output  power  of  the  projector. 


Figure  7.  An  example  photograph  of  the  oscilloscope  face  showing  the  direct 
arrival  and  the  (smaller)  reflected  arrival. 

Three  different  projectors  were  used;  a Dyna-Empire  model  J-9  was  used 
at  a spacing  of  100  m for  measurements  at  higher  acoustic  frequencies,  and  a 
more  powerful  Marine  Resources  J-ll  or  J-13  projector  was  used  at  a spacing  of 
300  m for  the  lower  frequency  measurements. 

The  reflection  coefficient  was  obtained  by  calculating  the  ratio  of 
the  amplitudes  of  the  reflected  pulse  and  the  direct  pulse.  Since  a ratio  of 
two  amplitudes  is  all  that  was  measured,  no  absolute  calibration  of  the 
equipment  was  necessary. 

The  difference  in  the  direct  and  reflected  path  lengths,  however, 
necessitated  the  use  of  a small  correction  to  compensate  for  the  different 
amounts  of  spherical  spreading.  Also,  at  higher  frequencies,  the  increased 
directivity  of  the  projectors  required  the  application  of  a correction  which. 


in  the  worst  case,  was  less  than  1 dB.  This  correction  was  made  by  assuming 
that  the  radiation  pattern  of  the  projector  was  that  of  a simple  circular 
piston.  The  far  field  (amplitude)  radiation  pattern  (normalized  to  1 at  0 = 0) 
for  a piston  is  given  by  (Urick,  1967) 


where  the  direction  6 is  measured  with  respect  to  the  normal  to  the  active 

surface  of  the  projector, 

a is  the  radius  of  piston  (2.85  cm  in  the  case  of  the  J-9 
and  5.1  cm  in  the  case  of  the  J-ll  and  J-13) 

X is  the  wavelength,  and 

is  the  first  order  Bessel  function  of  the  first  kind. 

The  angle  0 was  calculated  from  the  geometry  of  the  experiment  and  the  orien- 
tation of  the  projector.  The  J-9,  when  used,  was  oriented  to  point  horizon- 
tally in  the  direction  of  the  hydrophone;  the  J-ll  and  the  J-13,  on  the  other 
hand, were  suspended  to  point  vertically  upward. 

RESULTS 

Figure  8 through  11  show  the  corrected  reflectivity  in  decibels  as  a 
function  of  frequency.  Each  of  the  figures  shows  all  the  measurements  at  one 
location.  Each  diagram  indicates  the  horizontal  separation  and  depths  of  the 
projector  and  hydrophone,  the  grazing  angle  0,  and  the  type  of  projector.  One 
of  the  most  noticeable  features  of  the  plots  is  the  overall  low  reflectivity; 
values  of  -5dB  to  -lOdB  are  common.  The  other  remarkable  property  is  the 
large  variation  of  reflectivity  with  frequency;  sharp  minima  are  often  present. 
There  are  also  some  examples  showing  the  reflectivity  to  be  greater  than  0 dB 
suggesting  Interference  phenomena. 
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Figure  8.  Reflectivity  measured  in  Robeson  Channel  in  April,  1971 
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Figure  11.  Reflectivity  measured  in  Barrow  Strait  in  April,  1973. 
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Figure  11.  Continued. 
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Some  of  the  diagrams  occur  in  pairs  in  which  all  the  listed  parameters 
are  the  same.  The  difference  is  the  location  of  the  hydrophone  hole.  The 
two  holes  were  separated  by  an  azimuth  angle  of  about  20  degrees  in  order  to 
reflect  the  sound  from  different  samples  of  ice.  The  diagrams  which  are 
paired  in  this  way  are  Figures  9a  and  9b,  9c  and  9d,  9e  and  9f,  11a  and  lib, 
11c  and  lid,  and  lie  and  Ilf.  The  lack  of  detailed  agreement  between  the 
diagram  pairs,  expecially  those  in  Figure  11,  indicate  that  the  variations  in 
reflectivity  are  rot  mainly  due  to  the  gross  properties  of  the  ice  but,  rather, 
are  due  to  the  change  of  some  property  which  varies  rapidly  with  location. 

The  roughness  of  the  ice-water  interface  is  an  obvious  possibility. 

Another  interesting  feature  which  suggests  an  interference  phenomenon 
is  illustrated  in  the  photograph  shown  in  Figure  12. 


Figure  12.  This  photograph  shows  the  Interference  effect  at  the  beginning 
and  end  of  the  reflected  pulses.  This  effect  is  prominent  at 
frequencies  where  the  reflectivity  is  very  low. 
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This  picture,  which  was  taken  near  a reflectivity  minimum  shows  the  reflected 
pulse  sharply  decreasing  In  size  shortly  after  onset  and  then  Increasing  to  the 
original  amplitude  shortly  before  the  end  of  the  pulse.  As  will  be  shown, 
these  effects  are  generally  consistent  with  the  scattering  of  sound  from  a 
hard  Ice  surface  whose  RMS  roughness  is  only  a few  centimeters. 

SCATTERING  MODEL 

The  mathematical  model  used  to  establish  the  connection  between  our 
results  and  the  scattering  from  rough  surfaces  is,  in  essence,  an  application 
of  Kirchhoff's  Integral  Theorem  (Born  and  Wolf,  1959).  This  approach  to 
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where  k ■ — — with  f being  the  frequency  and  c the  speed  of  sound 

V2  2 

x +(dj+s)  the  distance  from  the  projector  to  the  element  of  area. 
2 2 

s ^ / (L-x  )+(d2+z)  the  distance  from  the  hydrophone  to  the  element 

‘ of  area. 

d^  ■ the  depth  of  the  projector 
d2  ■ the  depth  of  the  hydrophone 

z • the  height  of  the  element  of  area  above  the  mean  ice  surface 
L * the  horizontal  separation  of  the  projector  and  hydrophone. 

As  verification,  this  Integral  is  evaluated  in  the  Appendix  for  the 
case  of  a flat  surface  in  which  z(x)  “ 0 for  all  x.  A stationary  phase 
approximation  for  this  integral  yields 


♦(P)  * — 


lk(r  +s  ) 
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where  rQ  and  sq  are  the  values  of  r and  s corresponding  to  the  point  of 
geometric  reflection.  This  expression  is  what  one  would  expect  of  a 
perfectly  reflecting  hard  surface.  As  a further  check,  numerical  inte- 
grations yield  values  which  agree  well  both  in  magnitude  and  phase  with 
the  above  expression. 

The  remaining  problem  before  the  integral  can  be  evaluated  is  to 
simulate  the  rough  ice-water  boundary  with  a random  surface  z(x)  which 
has  a known  correlation  function  or,  equivalently,  a known  power  spectrum. 
To  calculate  this  spectrum,  a series  of  complex  Fourier  coefficients  of 
pre-asslgned  magnitudes  but  of  random  phases  was  generated  by  a computer 
and  then  the  inverse  Fast  Fourier  Transform  was  taken  to  provide  the 
random  surface  z(x).  The  mean-square  value  of  z(x)  is  easily  adjusted  and 
was  used  as  one  of  the  parameters  of  the  surface.  The  types  of  spatial 


power  spectra  that  were  tried  were:  "White"  with  a high  frequency  cut-off, 
l/f,  and  1/f  (where  f denotes  the  spatial  frequency  1/A).  Figure  13,  14 
and  15  show  examples  of  reflectivity  for  these  three  types  of  random  surface, 
the  reflection  geometry  being  the  same  and  the  RMS  value  of  each  surface  being 
2 cm.  It  can  be  seen  that  these  curves,  especially  Figure  14,  resemble  the 
experimental  curves  in  general  form.  This  agreement  strongly  suggests,  that 
the  observed  reflection  loss  was  primarily  the  result  of  scattering  from  the 
rough  surface. 

Using  the  ideas  developed  for  this  model,  it  is  possible  to  explain 
the  cancellation  phenomenon  that  was  illustrated  in  Figure  12.  Simply  ex- 
pressed, the  explanation  for  these  features  is  that  the  "Huygen"  wavelets  which 
arrive  first  at  the  hydrophone  come  from  the  area  of  minimum  travel  time,  the 
region  about  the  point  of  geometrical  reflection.  These  wavelets  produce  a 
large  signal  in  the  hydrophone.  Later  arriving  wavelets  fortuitously  have 
the  appropriate  amplitude  and  phase  to  cancel  this  almost  completely.  At  the 
end  of  the  pulse,  these  latter  wavelets  are  not  cancelled  by  the  initial  ones 
and  the  signal  becomes  large  again.  A calculation  of  the  reflected  pulse  was 
made  using  a formalism  similar  to  that  previously  described  except  that  the 
Green' 8 Function  was  now  time-dependent.  A calculated  plot  of  the  initial 
portion  of  the  pulse  is  shown  in  Figure  16.  As  can  be  seen,  the  amplitude 
decreases  as  the  contribution  from  a wider  and  wider  area  of  the  surface  is 
felt. 

An  aspect  of  the  problem  which  is  not  explained  by  this  model,  how- 
ever, is  that  the  observed  reflectivity  is  nearly  always  less  than  0 dB, 
whereas  the  calculated  reflectivities  are  often  positive.  This  implies  that 
some  of  the  energy  is  being  absorbed  via  a mechanism  not  Included  in  the 

model . This  could,  for  example,  be  the  absorption  of  sound  by  the  porous 
layer  at  the  water-ice  Interface. 
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Figure  15.  An  example  of  the  reflectivity  calculated  for  a random  surface 

whose  RMS  roughness  is  2 cm  and  whose  spatial  power  spectra  goes 
as  1/wavenumber2. 


CALCULATED  REFLECTION 


CONCLUSIONS 


The  acoustic  under-ice  reflectivity  measurements  obtained  during 
springtime  field  trips  in  the  years  1971,  1972  and  1973  are  presented.  The 
reflectivity  was  quite  variable  at  the  higher  frequencies  and  was  dependent 
on  the  detailed  structure  of  the  ice.  It  was  nearly  always  less  than  unity 
(0  dB)  for  all  angles  measured,  implying  a net  loss  of  energy  for  the 
reflection. 

In  a calculation  based  on  Kirchhoff's  Integral  theorem  the  ice 
surface  was  assumed  to  be  hard  and  the  random  undulations  of  the  ice-water 
interface  to  have  a RMS  value  of  2 cm.  This  calculation  predicted  variations 
in  the  reflectivity  which  agree  sufficiently  well  to  suggest  that  the  inter- 
ference mechanism  is  the  cause  of  the  variation.  The  model,  however,  does 
not  predict  a net  loss  of  energy;  reflectivity  predicted  by  the  model  is  often 
greater  than  unity. 
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APPENDIX  A 

This  appendix  derives  equations  I and  2 in  the  text. 

Assuming  single  frequency  acoustic  waves,  the  spatial  part  of  the 

wave  equation  reduces  to  the  Helmholtz  equation 

+ k^i|i  - 0 where  k ■ — 

c 

Kirchhoff's  Integral  Theorem  (Born  and  Wolf,  1959;  Baker  and  Copson  1953) 
states  that  for  a sufficiently  well  behaved  function  <J/,  the  value  of  at  a 
point  P can  be  expressed  as  an  integral  involving  the  function  and  its  normal 
derivative  evaluated  on  a surface  enclosing  the  point  P. 


ff 


elk8  3*  3 elks 

^ ds 


(A-l) 


where  s represents  the  distance  from  the  point  P to  the  surface  element  and 
0 

— represents  differentiation  with  respect  to  the  outward  normal. 
o n 

In  application  to  this  report,  the  enclosing  surface  consists  of 
the  ice-water  interface  and  a hemisphere  at  infinity.  An  additional  sphere 
enclosing  the  source  must  be  included  (and  integrate  over)  to  exclude  the 
discontinuity  at  that  point  from  the  enclosed  region.  The  integration  over 
this  small  sphere  yields  the  contribution  of  the  direct  arrival  and  will  be 


ignored  in  what  follows.  Also,  the  Integral  over  the  hemisphere  at  infinity 
is  assumed  to  be  zero.  Arguments  for  this  assumption  are  presented  in  Born 
and  Wolf  (1959).  This  leaves  only  the  Integral  over  the  ice-water  interface, 
and  this  will  yield  the  amplitude  of  the  reflected  wave. 


Figure  A-l:  A schematic  diagram  illustrating  the  notation  used  in  the 
derivation. 


With  reference  to  Figure  A-l,  let  r be  the  distance  from  the  source 


to  the  element  of  area  dS  and  s be  the  distance  from  this  element  to  the 


receiver  (point  P) . We  will  also  assume  that  the  undulations  of  the  ice  are 


gentle  enough  that  the  normal  derivative  is  approximately  equal  to  the 


derivative  with  respect  to  z.  In  the  absence  of  a surface,  the  value  of 

ikr 

0 

the  incident  wave  ip  at  the  element  of  integration  would  be  . The 


presence  of  a hard  surface  causes  the  value  of  the  function  to  double  and  the 


normal  derivative  to  be  equal  to  zero.  Equation  A-l  becomes 


ik(rfs)  d 


For  ease  of  computation,  let  us  now  assume  that  the  surface  irregu 


larltles  are  a function  of  x only,  that  is  z ■ z(x).  We  shall  Integrate  over 


y using  the  method  of  stationary  phase  (Born  and  Wolf,  1959,  page  749).  This 

•Lie  ( i } q \ 

uses  the  property  that  the  term  e varies  rapidly  with  y in  comparison 


to  the  rest  of  the  integrand  and  the  only  significant  contribution  to  the 


Integral  comes  from  the  vicinity  of  stationary  phase,  i.e 


(r+s)  - 0.  In  this  problem,  this  occurs  at  y ■ 0 


A statement  of  the  stationary  phase  approximation  is  given  by  the 


following  expression  wherein  the  point  of  stationary  phase 


ikh(y) 


In  the  present  derivation 


where 


We  have,  therefore,  using  equation  A-3  to  Integrate  over  y 


The  lengths  r and  s depend  on  the  surface  undulation  z which  must  be 


specified  as  a function  of  x in  order  that  the  integration  may  be  carried  out 


As  a check  on  the  equation  (A-4)  we  may  evaluate  the  Integral  for 


the  case  of  a flat  surface  (z  ■ 0).  The  method  of  stationary  phase  is  used 


again r where  this  time 
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